Density of States below the Thouless Gap in a Mesoscopic SNS Junction 
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Quasiclassical theory predicts an existence of a sharp energy gap Eg ~ hD/L^ in the excitation 
spectrum of a long diffusive superconductor-normal metal-superconductor (SNS) junction. We show 
that mesoscopic fluctuations remove the sharp edge of the spectrum, leading to a nonzero DoS for 
all energies. Physically, this effect originates from the quasi-localized states in the normal metal. 
Technically, we use an extension of Efetov's supermatrix cr-model for mixed NS systems. A non-zero 
DoS at energies E < Eg is provided by the instanton solution with broken supersymmetry. 



When a small piece of a normal (N) metal is placed in 
contact with a superconductor (S), paired electrons en- 
ter the N region, changing its excitation spectrum. How 
drastic are these changes? Recent studies demon- 
strate that the answer depends crucially on the type of 
dynamics in the N region: In the case of integrahle classi- 
cal dynamics, the density of states (DoS) of excitations is 
suppressed at low energies and vanishes nearly linearly at 
the Fermi level. Contrary, in N systems with chaotic dy- 
namics, coupling to a superconductor produces a gap in 
the DoS of the order of H/tc, where Tc is the typical time 
needed to establish contact with the superconductor. 

To illustrate this statement, consider a generic exam- 
ple of a chaotic NS system: an SNS junction made of a 
disordered conductor of size L connected to the S termi- 
nals. When the N metal is diffusive, with the mean free 
path I <^ L, and sufficiently long, with the Thouless en- 
ergy i?Th = hD/Lp' <C A (here D — vpl/S is the diffusion 
constant, and A is the superconductive gap in the ter- 
minals) , then Tc is given by the diffusion time across the 
N region, Tc ^ H/Etii. Thus, the energy gap Eg ^ E^h 
develops in the DoS ||,|. 

However, all the results mentioned above are based 
either on the quasiclassical theory of superconductivity 
and proximity effect ||-|^ or on the mead-field treatment 
of the random-matrix theory (RMT). Although usually 
this is a good approximation, it is interesting to check 
whether some effects, which are beyond quasiclassics, 
may lead to qualitative changes in the above picture. 

In this Letter we show that, indeed, mesoscopic fluctu- 
ations smear the hard gap in the quasiparticle spectrum 
of dirty SNS junctions, producing a tail of the subgap 
states with energies E < Eg. These low- lying states are 
due to the existence of quasi-localized states [| 10 in the 
N part of the junction which are weakly coupled to the S 
terminals, thereby having a larger Tc. The magnitude of 
this mesoscopic effect is controlled by the dimensionless 
normal-state conductance G = h/e^Rn of the junction, 
and is small (except for the vicinity of iJg) at G ^ 1. 
A related problem was considered recently in Rcf. pT] , 
where a hypothesis of universality and some nontrivial 
results from the RMT were used to find an expo- 
nentially small DoS below the mean-field gap in the N 
dot weakly connected to a superconductor. We will show 



that their conjecture holds for sufficiently narrow junc- 
tions close to Eg . Technically, our approach is completely 
different from Ref . as we employ the fully microscopic 
method based on the supermatrix tr-model method p^ ] 
for mixed NS systems |l^, which is free from any as- 
sumption about the random-matrix universality. We find 
a nontrivial instanton solution to the saddle-point equa- 
tions of the supermatrix cr-model and show that it is re- 
sponsible for the non-zero DoS at energies E < Eg. Our 
method is similar to the one used in |15|| to calculate the 
DoS in the quantum Hall regime. The relevance of the 
nontrivial saddle point solutions for the formation of the 
DoS tail in NS systems was pointed out in Ref. fl^ . 

We consider an SNS junction with the N region being 
a rectangular bar of size L x Ly x L^, coupled to the S 
terminals by the ideal contacts situated at a; = ±L/2. 
We neglect superconductivity suppression in the bulk 
terminals, provided that they are sufficiently large, and 
assume zero superconductive phase difference between 
them. The dimensionless conductance of the N region, 
G{L, Ly, Lz) = ATTvDLyLzl L, where v is the normal- 
metal DoS per a single projection of spin. 

Our results can be summarized as follows. The 
behavior of the subgap DoS {p{E)) depends on the 
relation between Ly,Lz, and the effective transverse 
length L±{E), with the latter scaling as L±{E 
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-(1 - S/Sg)(6-'i)/4, where Ly,,{E) = min(Lj,,„ Li(S)), 
and d is the effective dimensionality of the junction: 
d = for narrow junctions with Ly, Lz <C Lj_(E), d = 1 
for wider junctions with Ly ^ L±(E) ^ L^, and d = 2 
for films with Ly, Lz 3> L±(E). In the low-energy limit, 
E <^ Eg, the behavior of the DoS is log- normal in OD: 
lii{p{E)) ^ —Gln^{Eg/E), and is a power law in ID: 
In {p{E)) cx —ln{Eg/E). These results are similar to 
those for the distribution of relaxation times in ID and 
2D diffusive conductors although with an important 
difference: the basic energy scale in our case is given by 
Eg, whereas in Ref. it was the level spacing S = 1/vV . 

As a warm-up, we recall the standard quasiclassical 
approach to diffusive SNS junctions based on the 
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Usadel equation [|| for the quasiclassical retarded Green 
function G^{E, r). For the latter we will assume the an- 
gular parametrization, = cos + sin 9 cos </? + 
TySmOsmLp, where are the Pauli matrices in the 
Nambu space. In the absence of the phase difference be- 
tween the S terminals, one can set = 0, and the Usadel 
equation acquires the form (hereafter h = \) 



DV^0 + 2iEsme = 0, 



x = ±- 



L 



2' 



(1) 



The boundary conditions at the free edges of the N film 
are of von Neumann type. The local density of states is 
given by {piocaii^, E)) = 2z/Recos0(r). The substitution 
9{x) = 7r/2 -I- leads to the real equation for the 

function ip{x) with ip(±L/2) = 0. This equation can be 
easily integrated yielding the relation between E and the 
magnitude of '0(2;) at a; = 0: 
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V(o) 
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Th 



■\/sinh'0(O) — sinh-;/; 



(2) 



Eq. (H) has real solutions for ■0(0) only for E < Eg ~ 
ci?Th, with c = 3.12. At -E > Eg, V'(O) becomes complex 
resulting in the square-root singularity j^] in the DoS 
(here and below we provide results for the DoS integrated 
over the whole N region, {p{E)) = J dr (piocai(r, E)) ): 



= 2.30 S-\ E/Eg - 1 



(3) 



Below the mean field gap, at £; < Eg, Eq. (|) has 
two real solutions, '01(0) and V'2(0) > 0i(O), merging 
at E — Eg. They determine the corresponding solu- 
tions 6*1,2(2;) = '7r/2 + i0i,2(a;) to Eq. (|l|). Having real 
i/j{x), both of them do not contribute to the DoS at 
E < Eg. Usually the solution '0i(a;) is chosen by the 
continuity argument, as it obeys the natural condition 
lim^j^o 01(2;, E) = 0, while 0^2(2;, E) diverges in the limit 
of vanishing E. We will see however that it is this second 
solution which is responsible for the finite DoS below the 
quasiclassical gap. 

In order to extend the quasiclassical solution and take 
into account mesoscopic fluctuations, we will use the non- 
linear supermatrix a-model similar to that derived in 
. As we are interested in the average density of states 
only, it is sufficient to calculate the retarded single par- 
ticle Green function. The derivation of cr-model starts 
with representing in terms of the functional integral: 



G^(r,r',i?) 



bF{r)(l)^{r')e~^^'''^D(j)*D(l), (4) 



S[<j>]= J dr <j)+{r){E + iQ-'H)(f>{r). 

In this expression is the 4-component superfield consist- 
ing of commuting (complex) and anticommuting (Grass- 
mann) parts. The Hamiltonian 7i is a matrix in the 
Nambu-Gor'kov (N) space: 



r 



2m 



-Ef + U{r) 



(5) 



where A(r) = A6{\x\ — L/2), and U{r) is the random po- 
tential. In the absence of magnetic field, one has to dou- 
ble the field space in order to account the time-reversal 
symmetry, introducing the "time-reversal" (TR) space 
according to <I> = (0, iTy(j)*)'^ /^/2. This definition coin- 
cides with the one used in [16|; it differs from the nota- 
tions of by the factor iVy in the 0* sector. Pauli 
matrices operating in the TR space are denoted by <7i. 

After this definitions the derivation of cr-model is 
straightforward |13jl^Jl^. One has to carry out (i) av- 
eraging over the random potential with the correlator 
([/(r)C/(r')) = 6{r~r')/2T:i^T, (u) Hubbard-Stratonovich 
transformation introducing an 8 x 8 matrix Q acting in 
the product of FB, N, and TR spaces, (iii) expansion to 
the leading terms in VQ, E, and A. The result is 



(Aocai(r,S)) = JRc/ ^Qstr{fcAQ(r)}e 



-SIQ] 



(6) 



S[Q] 



-J dr str [^(Vg)^ -I- AiQ{iT^A + AE)] . (7) 



Here A ~ ctzTz, the matrix Q = U ^AU with the 
proper set of matrices U is subject to the condition 
Q = CQ^C^ where C = -T.:,(Tz[{l + k)a^ + {l-k)iay]/2, 
and k — diag(l, — 1)fb- The manifold of Q-matrices can 
be parametrized by 8 real and 8 Grassmann variables. 

The next step is to find the saddle point solutions to 
the action (|^) . We start from the simplest case when the 
Q matrix does not contain Grassmann variables. Then 
it splits into the Fermi-Fermi (FF) and Bose-Bose (BB) 
sectors which can be parametrized independently by four 
variables in each sector: 

Qff = Tz cos 0F [gz cos /cp + sin fcp {cFx cos xp + sin xf)] 
4- sin6'F(rj; cos</3f + Ty sinc^p), (8a) 

Qbb = [cz cos fee + Tz sinA;B(o'a; cosxb + (^y sIuxb)] 

X [tz cos 6b + ffz sin 9b[tx cos + Ty sin ^b)]- (8b) 

In terms of the angular parametrization (^, the action 
(^ acquires the form S = {ttv/2) / o?r(£FF — -Cbb), where 

£ff = D [(V6'f)^ + sin^ 6'f(V(^f)^ 

-I- cos^ 6'F(VfcF)^ -I- cos^ 6y sin^ fcF(VxF)^] 
+ AiE cos 6'f cos fcF — 4 A sin O-p cos (^f, 

Cbb = D [{VOBf + sin^ 0b(V</?b)' 
+ (Vfe)' + sin2 fcB(VxB)'] 
+ AiE cos 6b cos fee — 4 A sin cos fee cos ipB ■ 

The variables 0f,b and V5f,b coincide (at /cf.b = 0) with 
the standard Usadel angles, whereas their counterparts 
/cf,b and xf,b are new ingredients of the field theory. 
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Minimization of the action for a uniform superconduc- 
tor at <C A gives 9f,b = 7r/2 and A:f,b = 0, that 
provides the boundary conditions for Q in the N region. 
In the absence of the phase difference between the S ter- 
minals, one obtains v'f.b = fcp = and xf,b = const at 
the saddle point in the N part of the structure. Intro- 
ducing new variables ub = Ob + fee and /3b = — ^b in 
the BB sector, one obtains for the saddle-point action: 



5[0F,aB,/3B] - 2So[9f] - ^o[«b] - 5o[/3b], (9) 
T 



So[0] 



J dr [D{\/ef +4:iEcose] 



(10) 



Varying with respect to 6, one recovers Eq. (0) as the 
saddle-point equation for the action (|l0|). 

The Usadel equation (|l]) possesses, apart from the x- 
dependent solutions discussed above, solutions which de- 
pend on the transverse (y, z) coordinates. The role of 
the transverse dimensions will be discussed later, while 
now we will consider the OD case, relevant for sufRciently 
narrow junctions with Lx,Ly <C L±{E). Then, accord- 
ing to the previous analysis, the Usadel equation has 
two solutions, 9i{x) and 92{x). Therefore, the full ac- 
tion (||) has in total 8 different saddle point solutions: 
(6'f, ckB, /3b) = {0i,Oj,0k), with i, j, k = 1,2, that will 
be refered to as {i,j,k). However, only 4 of them with 
9-p = 9i can be reached by a proper deformation of the 
integration contour. 

The simplest is the supersymmetric saddle point 

(1.1.1) . In this case, Gaussian integrations over com- 
muting and anticommuting variables near it cancel each 
other, and the contribution to {p{E)) reduces to the form 
(^ with the vanishing DoS below Eg. Thus, the saddle- 
point approximation for the cr-model (|^) restricted to the 
supersymmetric saddle point is equivalent to the quasi- 
classical treatment based on the Usadel equation (|^). 

To get a non-zero DoS below Eg it is necessary to 
take saddle points with broken supersymmetry into ac- 
count. Such a solution with the lowest action is given by 

(1.1.2) [actually, a whole degenerate family of the sad- 
dle points, and, in particular, (1,2,1), can be obtained 
from it by rotation on the angle xb G [0, 27r)]. The 
key point is that Gaussian fluctuations near this sad- 
dle point have a negative eigenvalue which leads to an 
additional imaginary unity in the preexponent and, con- 
sequently, to the nonzero contribution to the DoS. This 
contribution is suppressed by the factor e~^'^, where 
AS* = — '5'o[6'2] > is the difference between 
the actions of the solutions ^1,2 (2;). Finally, the sad- 
dle point (1,2,2) has the action 2AS' and its contribu- 
tion can be disregarded at AS* ^ 1. Thus, the sub- 
gap DoS can be estimated with exponential accuracy as 



(piE)) 



-AS{E)_ j^qIqt^ calculate {p{E)) in 



the limiting cases E ^ Eg and E Eg. 

The solutions 6*1(2;) and d2{x) merge at Eg. This fact 
can be used to find the asymptotically exact result for 



{p{E)) in the energy range < 1 - E/Eg < 1. 

Let us start with the supersymmetric saddle point (1,1,1) 
and look at fluctuations around it. Almost all of them 
are hard, having a mass of the order of (or larger than) 
Eg. There are only 8 (corresponding to 4 commuting 
and 4 anticommuting variables) soft modes whose mass 
vanishes at E = Eg. Half of them transform the saddle 
point (1,1,1) to the instanton (1,1,2), and the other half 
transform it to the instanton (1,2,2). Below we will con- 
sider the case when AS 3> 1 that allows to disregard the 
contribution of the instanton (1,2,2) and to take Gaus- 
sian integrals over the corresponding soft and all hard 
fluctuations. As a result, we end up with 2 commut- 
ing (q and xb) and 2 Grassmann (^ and 5) variables 
parametrizing the relevant soft degrees of freedom in 
the matrix Q = e^"^^ /^e-^^ /^Ae"^" /^e^^ /\ Here the 
matrix W"" contains anticommuting variables: VFpB ~ 
(i/o(a;)/4)[(C + Oii-^y + -^zC^x) + (C - 0(*'ryCr^ + Tzi(Ty)\, 
W§p = Tx(Jx{Wy-qY' Txicfy, while in the absence of W"" , 
Q reduces to the form ^ with = /3b = ^1(2^), 
Qb — 0i{x) + iqfo{x), kp — ip-p = tpB = 0. The func- 
tion /o(x) is the normalized difference dip{x) = ip2{x) — 
■tjjiix) at Eg-. fo{x) = \imE^EgS'>p{x)/\\d'>p{x)\\, where 

||F(a;)||2 = (1/L) J^[^^^F'^{x)dx. Evaluating the action, 
integrating over the cyclic angle xb, and performing the 
x-integration, one obtains 



G 



S = G 
nc2E. 
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£ = 



2ci E^- E 



C2 



E, 



(11) 

(12) 



where c„ = j^^^^^cosh'4)[){x)fQ"' ^{x)dx/L, and?/jo(a;) — 
^pi^2{x. Eg); ci = 1.15, C2 — 0.88. The measure of in- 
tegration (including the Berezenian) is given by DQ — 
qdqdC,d£,. We take the contour of integration over q along 
the imaginary axis to get the convergent integral. Fi- 
nally, substituting / stT{kKQ)dV — 2iciV{4:,/e — q) into 
Eq. (^), we arrive at the onc-instanton correction to the 
quasiclassical result (||): 



ioo-t-O 



(piE)) 
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Im 



exp 



Gi^eq^ 



dq. (13) 



Eq. ( p^ is valid provided that e 3> G^"^!^ when the con- 
tribution of the instanton (1,2,2) can be neglected. Then 
the integral over q can be calculated by the saddle point 
method yielding at G"^/^ ^ e <C 1 in the OD case: 



{p{E 



Eg))oi) - 4^ 



G^ 



exp 



iG£3/2). (14) 



This result can be generalized for a normal dot of an 
arbitrary shape coupled to a superconductor (cf. |p],pT[), 
provided that the numbers c, Cn are defined with the use 
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of the exact solutions 0i.2(r) of the Usadel equation in 
a given geometry: c„ = (l/V) J cosh V;o(r)/Q"~^(r) dr. 
The functional form of the result ( p^ coincides with the 
RMT conjecture of Ref. An exact correspondence 

with the RMT prediction is expected in the case of 
weakly transparent NS interfaces. 

Coming back to the planar SNS junction, we turn 
to the limit of small energies, E <^ Eg. Here one has 
i'li.x) ~ and ip2{x) ~ A{1 — 2\x\/L), where, according 
to Eq. (|), A « \nlEg/E) (cf. Ref. [||). Thus, the action 
of the instanton (1,1,2) becomes AS{E) w —So[92] = 
TrEThA"^ /6, and the result for DoS reads 



{p{E « Eg)) 



OD 



■ cxp 



-^In^:^ 

4 E 



(15) 



Now let us consider the role of the saddle-point solu- 
tions which depend on the transverse coordinates y, z. 
At E Eg, one has to retain only soft modes associated 
with the instanton (1,1,2). As a result, the action ( |ll| ) 
acquires a gradient term: 



S^G I ^^IJ^(^V^q)' + V-eq'-^ 



(16) 



where the Grassmann variables are discarded as we are 
not interested in the preexponent. Comparing the first 
and the second term in Eq. (p^), one extracts the char- 
acteristic transverse scale L±{E) = L{cc2)^^^^e^^^^ ~ 
Le~^/'*, which determines the effective dimensionality of 
the system. If Ly or Lz is shorter than L±{E), then 
it costs too much energy to have gradients in that direc- 
tion, and the corresponding dimension "freezes out" . The 
OD case considered above referred to the limit Ly, Lz ^ 
L±{E). Otherwise, an instanton will appear in the trans- 
verse direction to minimize the total action. In the ID 
case {Ly ^ L±{E) ^ Lz), the action ( |l^ ) achieves its 
stationary point at q{y) ~ 3a/£ cosh~^(y/L^), leading to 



{p{E ^ Eg)) 



SI/I'D 



1 / 127r , 

- cxp —Jcc2vDLz e ' 

V 5 



(17) 



Analogously, in the 2D case {Ly,Lz ^ L^[E)), the 
quasiparticle DoS tail has the form 



{p{E ^ Eg)) 



1 



2D 



exp(-24.4j/DLe) 



The length L±{E) diverges at Eg, indicating that any 
junction becomes effectively OD close to the quasiclassi- 
cal gap. However, different parts of the DoS tail may 
exhibit different exponents, from 1 to 3/2 [cf. Eqs. (14), 
(|l|), and (|l|)], with some cross-over in between. 

The value of L±{E) is getting shorter as E decreases, 
and becomes of the order oi L a.t E < Eg. The solution 
of the ID problem a.t E <^ Eg can be found following 
Ref. 1^. The function 4'2{x,y) has a sharp peak at the 
center of the instanton and with the logarithmic accuracy 



is given by ip2{x,y) = 
for the DoS then reads 



-41n(2v/a;2 + J/7L). The resuh 



{p{E « Eg)) 



ID 



(19) 



The DoS a,t E <^ Eg in the 2D case can be calculated in 
the same way, and appears to be vanishingly small. 

To conclude, we have shown that mesoscopic fluctua- 
tions smear the quasiclassical gap in the DoS of a diffusive 
SNS junction. The tail in the DoS is due to the states 
anomalously localized in the N part of the junction and 
weakly coupled to the S terminals. Technically, the tail 
is described by instantons with broken supersymmetry. 

At the final stage of the preparation of this manuscript 
we became aware of Ref. jT^ where a similar instanton 
solution is found for the problem of tail states below Eg 
in a superconductor with magnetic impurities. 
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